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Abstract 

Taste receptor genes are functionally important in animals, with a surprising exception in the bottlenose dolphin, which shows 
extensive losses of sweet, umami, and bitter taste receptor genes. To examine the generality of taste gene loss, we examined seven 
toothed whales and five baleen whales and sequenced the complete repertoire of three sweet/umami (T1Rs) and ten bitter (T2Rs) 
taste receptor genes. We found all amplified T1Rs and T2Rs to be pseudogenes in all 1 2 whales, with a shared premature stop codon 
in 1 0 of the 1 3 genes, which demonstrated massive losses of taste receptor genes in the common ancestor of whales. Furthermore, 
we analyzed three genome sequences from two toothed whales and one baleen whale and found that the sour taste marker gene 
Pkd2H is a pseudogene, whereas the candidate salty taste receptor genes are intact and putatively functional. Additionally, we 
examined three genes that are responsible for taste signal transduction and found the relaxation of functional constraints on taste 
signaling pathways along the ancestral branch leading to whales. Together, our results strongly suggest extensive losses of sweet, 
umami, bitter, and sour tastes in whales, and the relaxation of taste function most likely arose in the common ancestor of whales 
between 36 and 53 Ma. Therefore, whales represent the first animal group to lack four of five primary tastes, probably driven by the 
marine environment with high concentration of sodium, the feeding behavior of swallowing prey whole, and the dietary switch from 
plants to meat in the whale ancestor. 
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Introduction 

Taste perception is fundamental for the survival of animals 
ranging from insects to mammals (Yarmolinsky et al. 2009). 
The sense of taste is highly specialized to sense nutritious or 
harmful compounds in potential food sources and is therefore 
essential to trigger or regulate feeding behaviors in animals 
(Bachmanov and Beauchamp 2007). Animals are commonly 
believed to have five primary taste modalities: Sweet, umami, 
bitter, sour, and salty in mammals; and sweet, bitter, water, 
carbonation, and salty in insects (Bachmanov and Beauchamp 
2007; Yarmolinsky et al. 2009). Although mammals and in- 
sects diverged from their common ancestor approximately 
900 Ma (Hedges et al. 2004), both groups of animals share 
the same fundamental principles for encoding taste informa- 
tion (Yarmolinsky et al. 2009). 



Taste is activated by the physical interactions between 
tastant molecules from the external environment and taste 
receptors, the latter of which have been decoded in the last 
decade (Bachmanov and Beauchamp 2007; Yarmolinsky et al. 
2009). Vertebrate taste receptors are either ion channels (sour 
and salty) or two families of G protein-coupled receptors (T1 Rs 
for sweet/umami and T2Rs for bitter), whereas gustatory re- 
ceptors (GRs) are primarily responsible for taste in insects. 
Because of the essential roles of taste, taste receptor genes 
are believed to be indispensable from insects to mammals. For 
example, the sweet/umami taste receptor gene repertoire 
(T1Rs) generally contains three members across mammals; 
the number of GRs remains similar in most insects (Nei et al. 
2008). Although vertebrate bitter taste receptor genes (T2Rs) 
vary remarkably in number from 3 in chicken to 49 in frog (Shi 
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and Zhang 2006), the bitter taste is still required to detect 
toxins in food sources for these animals. However, numerous 
pseudogenizations of taste receptor genes have been discov- 
ered in animals along with an increasing number of available 
genome sequences. For example, the sweet taste receptor 
gene (T1R2) is lost in the chicken, zebra finch, cat, vampire 
bats, and western clawed frog (Li et al. 2005; Shi and Zhang 
2006; Zhao et al. 2010); the umami taste receptor gene T1R1 
is lost in the giant panda, western clawed frog, and bats (Shi 
and Zhang 2006; Zhao, Yang, et al. 2010; Zhao et al. 2012). 
Strikingly, the draft dolphin genome (2.59x coverage) lacks 
intact genes responsible for sweet, umami, and bitter tastes 
(Jiang et al. 2012). The lack of bitter taste is particularly unex- 
pected because natural toxins typically taste bitter, and bitter 
taste thus represents an important natural defense against the 
ingestion of poisonous chemicals from the external environ- 
ment (Glendinning 1994). 

In addition to taste receptors, signaling pathways down- 
stream of taste receptors are also essential for taste function. 
Specifically, three ion channels (TRPM5, PLCP2, and CALHM1) 
have been identified as key components of taste signal trans- 
duction to detect sweet, umami, and bitter tastes (Zhang et al. 
2003; Taruno et al. 2013). TRPM5 (Transient receptor poten- 
tial cation channel subfamily M member 5), a taste-specific 
TRP ion channel, is coexpressed with T1 R and T2R cells in taste 
buds, and mice deficient for Trpm5 exhibit abolished or se- 
verely reduced sensitivity to sweet, umami, and bitter stimuli 
(Zhang et al. 2003; Damak et al. 2006). Similar to TRPM5, 
PLCP2 (phospholipase Cp2) shares overlapping expression 
patterns with TRPM5, and its knockout mice show a selective 
and complete loss of sweet, umami, and bitter responses 
(Zhang et al. 2003). CALHM1 (calcium homeostasis modulator 
1) contributes to the neurotransmission of taste stimuli; the 
loss of CALHM1 rendered severely impaired responses to 
sweet, umami, and bitter tastants (Taruno et al. 2013). 
Thus, Trpm5, Plc(52, and Calhml likely undergo the relaxation 
of selective constraints if sweet, umami, and bitter tastes have 
been lost for a considerable amount of time. 

To examine the generality of taste loss in whales and de- 
termine when whales lost tastes, we sequenced T1Rs and 
72/?s in 11 species of whales: The bottlenose dolphin 
{Tursiops truncates), pilot whale (Globicephala melas), white- 
beaked dolphin {Lagenorhynchus albirostris), Atlantic white- 
sided dolphin (L. acutus), finless porpoise {Neophocaena 
phocaenoides), sperm whale (Physeter catodon), Bryde's 
whale (Balaenoptera edeni), Omura's whale (B. omurai), fin 
whale (B. physalus), minke whale (B. acutorostrata), and bow- 
head whale {Balaena mysticetus) (fig. 1). Our sample con- 
tained representatives from both major lineages of whales: 
Odontoceti (toothed whales) and Mysticeti (baleen whales) 
(McGowen et al. 2009) (fig. 1). We included three sweet/ 
umami taste receptor genes (T1R1-3) and all ten bitter taste 
receptor genes (T2R1 -3, T2R5, T2R16, T2R38-39, T2R60, and 
T2R62a-62 b) identified from the draft dolphin genome (Jiang 



et al. 2012). We also examined the draft genome sequences 
of the bottlenose dolphin, Yangtze River dolphin (Zhou et al. 
2013), and minke whale (Yim et al. 2014) to identify genes 
potentially responsible for sour and salty tastes. To test 
whether taste signaling pathways have been relaxed from 
selective constraints because of taste loss, we analyzed the 
Calhml, Plc/32, and Trpm5 sequences from several whales in 
comparison with other mammals. We show that the sweet, 
umami, bitter, and sour taste receptor genes are pseudogen- 
ized, whereas salty taste receptor genes are evolutionarily con- 
served in toothed and baleen whales. We further show the 
relaxation of selective constraints on taste signaling pathways 
along the ancestral branch leading to whales (P/c/32 and 
Calhml) or after the divergence of whales (TrpmS). 

Materials and Methods 

Polymerase Chain Reaction and DNA Sequencing 

By searching the draft genome sequence of the bottlenose 
dolphin (2.59X), we identified all three members of T1Rs, all 
ten members of T2Rs, and Calhml. These sequences were 
used to design a suite of primers (supplementary table S1, 
Supplementary Material online), which were employed to am- 
plify the 14 genes in other whales. Genetic material was ob- 
tained from the Institute of Hydrobiology, Chinese Academy 
of Sciences, and the marine mammal tissue archive held at 
Zoological Society of London. Our samples contained six spe- 
cies of toothed whales (bottlenose dolphin, pilot whale, white- 
beaked dolphin, Atlantic white-sided dolphin, finless porpoise, 
and sperm whale) and five species of baleen whales (Bryde's 
whale, Omura's whale, fin whale, minke whale, and bowhead 
whale) (fig. 1). The genomic DNAs were isolated using the 
Qiagen DNeasy kit. The polymerase chain reaction (PCR) mix- 
tures (30 jil) contained 0.5 pj genomic DNA (1 0 ng/jil), 1 5 \x\ of 
2x PCR solution (Takara Premix Taq), and 1.2pJ of each 
primer (10jiM). When the above PCRs did not work, we 
used the following PCR mixtures (25 jil): 0.5 jil genomic 
DNA (10ng/jil), 2.5 jil of 10x PCR buffer, 2.5 jil dNTPs, 
1.5 jil MgS0 4 , 0.75 jil of each primer (10 jiM), 0.5 jil (1U/ jil) 
of KOD-Plus-Neo DNA polymerase (Toyobo), and 16 jil H 2 0. 
All PCRs were conducted on a BioRadTI 00 Thermal Cycler 
under the touchdown conditions as follows: 5min of initial 
denaturation, 15 cycles of denaturation at 95 °C for 30 s, 
65 °C for 30s (the temperature was decreased by 1 °C per 
cycle), extension at 72 °C for 60s; followed by 25 cycles with 
95 °C for 30s, 50 °C for 30s, 72 °C for 60s; and a final ex- 
tension at 72 °C for 5min. Sequencing reactions were per- 
formed directly from both strands with the same primer sets 
as those used for the PCR amplifications. When the direct 
sequencing did not work, the PCR products were cloned 
into the pMD19-T vector (Takara) and transformed to DH- 
5a competent cells. Three to five positive clones of each 
PCR product were screened and sequenced. All newly 
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Fig. 1. — The species tree of whales with common mutations of T1Rs and T2Rs. The tree topology and divergence dates follow previous studies 
(Hasegawa etal. 2003; McGowen etal. 2009). Branch lengths are not drawn to scale; the numbers at the nodes are the divergence times in millions of years. 
T1R1 and T1R3 encode the umami taste receptor, T1R2 and T1R3 encode sweet taste receptors, and T2Rs confer bitter taste. For each gene, shared frame- 
shifting indels or premature stop codons across toothed and baleen whales are indicated by "S," shared frame-shifting mutations in one of the two major 
lineages of whales are indicated with "S*" (toothed whales) or "S # " (baleen whales), and unshared ones are shown with "U." Additionally, "+" indicates 
the presence of an intact ORF, " - " denotes the absence from the genome, while no signs in any gene indicate no amplifications despite numerous attempts. 



generated sequences by PCRs were submitted to GenBank 
under accession numbers KJ547495-KJ547591 . 



Data Mining from Mammalian Genome Sequences 

For the identification of multiple-exon genes, we used mouse 
genes as query sequences to conduct TBLASTN searches 
(Altschul et al. 1990) to identify sour taste marker genes 
{Pkd2l1 and Pkd1l3), the candidate salty taste receptor genes 
(Scnnla, Scnnlb, and Scnnlg), and three taste signaling path- 
way genes (Calhml, Trpm5, and Plc/32) from three publicly 
available genome sequences of whales (bottlenose dolphin, 
Yangtze River dolphin, and mink whale). To validate the 
BLAST results, we retrieved the matched scaffolds from 
genome sequences, and the exon/intron structures were de- 
termined using a combination of GeneWise (Birney et al. 
2004) and NCBI-BLAST2 (Altschul et al. 1997) programs. For 
the identification of single-exon genes (T2Rs), we used all T2Rs 
from human, rat, dog, and chicken as queries to TBLASTN the 
three whale genomes following a previous study (Shi and 
Zhang 2006) and confirmed the presence of seven 



transmembrane domains using the TMHMM method 
(Sonnhammer et al. 1998). All candidate genes were verified 
by the best hits with known genes of interest using BLASTP 
searches against the entire GenBank (Shi and Zhang 2006). All 
sequences culled from the three whale genomes are provided 
in the supplementary data set S1, Supplementary Material 
online. 

Evolutionary Analysis 

Nucleotide and protein sequences were aligned with 
CLUSTALX (Thompson et al. 1997) and modified by eye 
with Bio-Edit (Hall 1999). The best-fit substitution models 
were selected using the jModelTest2 program for each data 
set according to the Akaike information criterion and the 
Bayesian information criterion (Darriba et al. 2012). 
Maximum likelihood (ML) and Bayesian phylogenetic trees 
were reconstructed with PhyML version 3.1 (Guindon et al. 
2010) and MrBayes version 3.1.2 (Huelsenbeck and Ronquist 
2001), respectively, following a previous study (Lin et al. 
2014). PhyML 3.1 used the nearest neighbor interchange 
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algorithm to search ML trees with 100 bootstrap replicates, 
whereas MrBayes 3.1 .2 ran 1 million generations to reach the 
convergence with the standard deviation (SD) of split frequen- 
cies lower than 0.01 . The ancestral sequences of whales were 
reconstructed using the Bayesian approach implemented in 
the BASEML program from the PAML package (Yang 2007). 
Nonsynonymous (c/ N ) and synonymous (d s ) distances and their 
standard errors from between sequences were estimated 
using the modified Nei-Gojobori method (Zhang et al. 
1998) implemented in MEGA5 (Tamura et al. 2011). The 
ML estimates of the ratios of nonsynonymous and synony- 
mous nucleotide substitution rates ((£> = d N /d s ) were gener- 
ated by the CODEML program in PAML (Yang 2007). 

Results 

Evolution of Taste Receptor Genes in Whales 

To determine how widespread the taste loss is in cetaceans, 
we examined all taste receptor genes that are responsible for 
five primary taste modalities in mammals: Sweet, umami, 
bitter, sour, and salty. First, we sequenced all members of 
the sweet/umami (3 T1Rs) and bitter (10 T2Rs) receptor 
gene families in six toothed whales and five baleen whales, 
and we also obtained these genes in a fully sequenced 
genome of one additional toothed whale {Lipotes vexillifer). 
Second, we analyzed three whole genome sequences of two 
toothed whales (7". truncates and Li. vexillifer) and one baleen 
whale (B. acutorostrata) to survey the sour taste marker gene 
(Pkd2l1) and the candidate salty taste receptor genes {Scnnla, 
Scnnlb, and Scnnlg). 

Pseudogenization of Sweet/Umami (T1Rs) and Bitter 
(T2Rs) Taste Receptor Genes 

Mammals typically have only three members of T1R genes 
(T1R1, T1R2, and T1R3), among which T1R1 and T1R3 
encode two subunits of the umami taste receptor, whereas 
T1R2 and T1R3 encode two proteins that function as the 
sweet taste receptor (Nelson et al. 2001, 2002; Sainz et al. 
2001 ; Li et al. 2002; Zhao et al. 2003). T1Rs include six exons: 
The first five encode the extracellular domain, whereas the 
sixth exon encodes all remaining domains of the receptors 
(Li et al. 2002). We successfully amplified T1R1 (exon 6), 
T1R2 (exons 3 and 6), and T1R3 (exon 3) from six toothed 
whales and five baleen whales. All newly obtained sequences 
were characterized by frame-shifting indels (deletions or inser- 
tions) and/or premature stop codons, which are hallmarks of 
pseudogenes (fig. 1 and supplementary table S2, 
Supplementary Material online). This finding suggests that 
all T1Rs in the examined whales are pseudogenes. 
Moreover, we identified one 1 -bp deletion of T1R3 and one 
premature stop codon of T1R1 that are shared between 
toothed and baleen whales, which are indicative of an ances- 
tral pseudogenization in the common ancestor of whales. 



However, we failed to identify any common disruptive muta- 
tions from both exons 3 and 6 of T1R2 (fig. 1 and supplemen- 
tary table S2, Supplementary Material online). Because T1R3 
encodes the shared subunit of the sweet and umami taste 
receptors, our results strongly suggest that sweet and 
umami tastes were lost in the common ancestor of whales. 

Bitter taste receptors are encoded by T2Rs, which lack 
introns and are -900 bp in length (Adler et al. 2000; 
Chandrashekar et al. 2000). We attempted to sequence 
ten T2Rs that are orthologous to those identified from 
the dolphin genome (Jiang et al. 2012). With one excep- 
tion, all amplified T2Rs have disrupted open reading frames 
(ORFs) characterized by frame-shifting indels and/or prema- 
ture stop codons (fig. 1 and supplementary table S2, 
Supplementary Material online). The exception is T2R16, 
which is intact in the three baleen whales (fig. 1). We 
also computed the nonsynonymous (c/ N ) to synonymous 
(d s ) distances for T2R16 using MEGA5 (Tamura et al. 
2011); the c/ N /c/ s ratio does not significantly differ from 1 
in all pairwise comparisons of the three species (mean±SD, 
0.48 ±0.07; P>0.1, Z test). In support of this finding, we 
estimated the mean d N /d s ratio of T2R16 pseudogenes 
(mean±SD, 0.70 ±0.35) and that of other T2R pseudo- 
genes in baleen whales (mean±SD, 1.37 ±0.78), both 
mean c/ N /c/ s ratios are not significantly different from 
0.48, the mean d^d s ratio of T2R16 intact genes 
(P>0.05, two-tailed f-test). Thus, T2R16 appears to un- 
dergo relaxation of functional constraints in the three 
baleen whales, which is comparable to the pseudogeniza- 
tions in two additional baleen whales and five toothed 
whales (fig. 1). Among the remaining nine T2R pseudo- 
genes, eight contain shared frame-shifting indels and/or 
premature stop codons between toothed and baleen 
whales, whereas one (T2R5) includes these pseudogene 
hallmarks shared in toothed whales only. The pseudogen- 
ization of T2R5 may have occurred in baleen whales earlier 
than toothed whales because we failed to obtain this gene 
from PCR amplifications for five baleen whales and from 
the mink whale genome (fig. 1 and supplementary fig. 
S1A Supplementary Material online). Indeed, although 
T2R5 is flanked with PRSS37 and T2R3 in the bottlenose 
dolphin and Yangtze River dolphin genomes, we found 
that the two genes are adjacent to each other in scaf- 
fold_20 of the minke whale genome, suggesting a true 
loss of T2R5 in the minke whale. In addition, we identified 
comparable numbers of T2Rs in the three whale genomes 
(10 in bottlenose dolphin, 11 in Yangtze River dolphin, and 
13 in mink whale), and we were able to find all ten T2Rs 
in the Yangtze River dolphin and nine of ten T2Rs in the 
mink whale (supplementary fig. S1/\, Supplementary 
Material online), suggesting that the ten T2Rs could repre- 
sent the bitter taste receptor gene repertoire in all whales. 
All T2Rs identified from the three whale genomes are 
pseudogenes, except the intact T2R16 in the minke 
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whale. Together, these results strongly indicate that bitter 
taste has been lost in all whales examined, and the relax- 
ation of selective constraint on bitter taste receptor genes 
most likely occurred in the common ancestor of whales. 

Given that toothed and baleen whales diverged 36 Ma 
and whales and even-toed ungulates diverged 53 Ma 
(Hasegawa et al. 2003; McGowen et al. 2009), our genetic 
evidence suggests that the pseudogenizations of sweet, 
umami, and bitter taste receptors took place in the 
common ancestor of whales between 36 and 53 Ma 
(fig. 1), probably shortly after 53 Ma, because we observed 
a number of common frame-shifting mutations in many 
genes predating the divergence of both major lineages of 
whales (fig. 1 and supplementary table S2, Supplementary 
Material online). 

Pseudogenization of Sour Taste Marker Gene (Pkd2l1) 

Two sour taste marker genes (Pkd1l3 and Pkd2l1), the ex- 
pression of which is indispensable in sour taste functioning, 
encode the transient receptor potential channel members 
that were proposed to function as candidate sour taste 
receptors (Huang et al. 2006; Ishimaru et al. 2006). We 
used mouse genes as a query to search against the draft 
genome sequences of the dolphin, Yangtze River dolphin, 
and minke whale (Zhou et al. 2013; Yim et al. 2014). We 
failed to identify unambiguous sequences of Pkd1l3 due to 
incomplete sequencing but instead obtained nearly com- 
plete coding sequences of Pkd2l1 (supplementary table 
S2, Supplementary Material online). After aligning the 
three newly identified sequences with mouse Pkd2l1 
(GenBank accession no. NM_181422), we found multiple 



premature stop codons in each of the three whales with 
available genome sequences (fig. 2 and supplementary 
table S2, Supplementary Material online), which suggested 
that none of these sequences are functional. In the dolphin 
and minke whale sequences, the first premature stop 
codon is located at exon 4, leading to the loss of func- 
tional domains from exon 5 to exon 15. In the sequence of 
the Yangtze River dolphin, the 5 r -most premature stop 
codon resides at exon 5, rendering the remaining ten 
exons nonfunctional. Furthermore, the dolphin and minke 
whale share the first premature stop codon located at exon 
4, which indicates the possible pseudogenization of Pkd2l1 
in the common ancestor of whales, although exon 4 of the 
Yangtze River dolphin is complete and intact (fig. 2 and 
supplementary table S2, Supplementary Material online). In 
addition, two common premature stop codons were found 
in the dolphin and Yangtze River dolphin (fig. 2 and sup- 
plementary table S2, Supplementary Material online), sug- 
gesting a loss of function for Pkd2l1 that predates the 
divergence of the two whales after separation from the 
minke whale. To determine when the functional relaxation 
of Pkd2l1 occurred, we also inferred the Pkd2H sequence 
of the common ancestor of the three whales (blue circle in 
supplementary fig. SIB, Supplementary Material online) 
using a Bayesian approach (Yang et al. 1995) and esti- 
mated the co ratios of Pkd2l1 for the common ancestor 
and 16 nonwhale mammals (supplementary fig. S1£, 
Supplementary Material online). We found that a model 
that allows a variation in co between the common ancestor 
(co 1= 0.15) and all 16 nonwhale mammals (co 2 = 0.35) is 
significantly better than a simpler model that assumes the 
same co (co 0 = 0.16) across the tree (P=0.016), which 



Pkd2ll Exon 4 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(597) 5' ttgcggcagctgcgcgtgcgcaatgactcctgtgtggttcatgaagacttccgggaggac 3' 

(597) 5' ctgtggcagctgagggtccgcaatgactccBtggtgcacgaggacttccatgaggac 3' 

(597) 5' ctgtggcagctgagggtctgcaatgactccggagtggtgcacgaggacttccacgaggac 3' 

(597) 5' ctgcagcagctgagggtccgcaatgactccfgabtggtgcacgaggacttctgcgaggac 3' 



PkdUl Exon 5 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(733) 5' TGGACATACCATTCCCAGAATGAGCTGGGTGGCTCCTCCCACTGGGGCAGGCTCACAAGC 3' 

(733) 5' gACCTACCACHCAGGATGAGCTGGGGGGCTCCTCTCACTGGGGCCAGCTCACAAGC 3' 

(733) 5' gACCTACCACTCGCAGGATGAGCTGGGGGGCTCCTCTCACTGGGGCCAGCTCACAAGC 3' 

(733) 5' TGGACCTACCACTCGCAGGATGAGCTGGGGGGCTCCTCTCACTGGGGCCGGCTCACAAGC 3' 



PM2U Exon 12 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(1900) 5' GAGATCACCGCTGCCTTCACCAGGTTTGATCAGGATGGGGACCACATACTGGATGAGGAG 3' 

(1900) 5' GAGCTCACGGCTGCCTTCACCAGGTTTGACgGATGGGAATCACATCCTGGATGAGAAA 3' 

(1900) 5' GAGCTCACGGCTGCCTTCACCAGGTTTGA(@GATAGGAATCACATCCTGGATGAGAAA 3' 

(1900) 5' GAGCTCACGGCTGCCTTCACCAGGTTTGGCCAAGATGGGACTCACATCCTGGATGAGAAG 3' 



Fig. 2. — An alignment of Pkd2l1 containing the ORF-disrupting mutations in two toothed whales (bottlenose dolphin and Yangtze River dolphin), one 
baleen whale (mink whale), and the outgroup (mouse). Codons in the correct reading frame are indicated by shading, and premature stop codons are boxed; 
the numbers in parentheses indicate the nucleotide positions following the mouse sequence. 
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suggested that a functional relaxation of Pkd2l1 arose in 
the common ancestor of whales. 

Purifying Selection on Candidate Salty Taste Receptor 
Genes (Scnn1a f Sen n 1b, and Scnnlg) 

The epithelial sodium channel ENaC has been proposed to 
function as a candidate salty taste receptor, which usually in- 
volves three subunits encoded by Scnnla, Scnnlb, and 
Scnnlg (Heck et al. 1984; Canessa et al. 1994; Firsov et al. 
1998; Yarmolinsky et al. 2009). We examined the draft 
genome sequences of the dolphin, Yangtze River dolphin, 
and minke whale. In contrast to the widespread pseudogen- 
ization of other taste receptor genes, we were able to identify 
three salty taste receptor genes with nearly complete coding 
regions and intact ORFs (supplementary table S2, 
Supplementary Material online), which suggested that salty 
taste is retained in the three whales. Consistent with this find- 
ing, the c/ N /c/ s ratio of each of the three genes was significantly 
lower than 1 in all pairwise comparisons of the three species 
(P< 0.001, Ztest). These results indicate that the candidate 
salty taste receptor genes have experienced strong purifying 
selection, and salty taste may be functional in whales. 

Relaxation of Selective Constraint on Taste Signaling 
Pathways in Whales 

To test whether widespread losses of taste receptor genes 
have affected the evolution of taste signaling pathways, we 
examined three genes {Calhml, Trpm5, and Plcfi2) that are 
known to be involved in taste signaling (Zhang et al. 2003; 
Taruno et al. 2013). Specifically, Calhml encodes the 
CALHM1 ion channel that is required for the neurotransmis- 
sion of sweet, bitter, and umami tastes (Taruno et al. 2013); 
Trpm5 and Plc/32 encode a TRP ion channel and a 
phospholipaseC, respectively, both of which are common sig- 
naling molecules for sweet, umami, and bitter taste transduc- 
tion (Zhang et al. 2003). 

We sequenced all two coding exons of Calhml from six 
toothed whales and four baleen whales and identified an ad- 
ditional Calhml from the genome of one toothed whale 
(Yangtze River dolphin) (figs. 3 and 4). With two exceptions, 
the ORF of Calhml is intact in nine other sequenced whales, 
suggesting a functional role in most whales. The first excep- 
tion is one toothed whale (Delphinapterus leucas), which 
showed one 1 -bp deletion in exon 1, whereas the other ex- 
ception is a baleen whale (B. omurai), where one 1 -bp inser- 
tion in exon 2 occurred (fig. 3). These indels result in a shift in 
the ORF and could create multiple premature stop codons, 
which are indicative of a loss of function. Moreover, the lack 
of common ORF-disrupting indels in the two whales suggests 
independent events of the relaxation of functional constraint. 
To determine whether the relaxation of functional constraint 
occurred before or after the radiation of extant whales, we 
estimated the co ratios of Calhml in whales using a likelihood 



approach. We computed likelihood tests on two data sets. 
Data set I included 26 nonwhale mammals and one ancestral 
sequence of all whales, whereas data set II contained all mam- 
mals (26 nonwhale mammals and 1 1 whales) after the re- 
moval of premature stop codons (table 1). First, we 
examined data set I under the assumption of a uniform co 
across the tree (model A in table 1); co was estimated to be 
0.06, which suggested an overall purifying selection on 
Calhml in mammals. Assuming that the ancestral sequence 
of all whales has co 2 and other branches have (model B in 
table 1), we estimated co 2 = 0.17, and model B fits the data 
significantly better than model A (P= 0.021), indicating the 
common ancestor of all whales has a significantly higher co 
than other mammals. This result supports the hypothesis stat- 
ing that a relaxation of functional constraint started from the 
common ancestor of all whales. Indeed, the functional relax- 
ation of Calhml is well reflected by the elevated branch length 
of the ancestral branch leading to whales (fig. 4). Second, we 
analyzed data set II and allowed a variation in co between the 
stem whale branch and all branches connecting 1 1 whales 
(model D in table 1). The stem whale branch was estimated 
to have a significantly lower co than the other whale branches 
(P=0.01) after comparing model D with model C, the latter 
of which is a simpler model assuming the same co for all whale 
branches (table 1 ). This analysis suggests the further relaxation 
of functional constraint after the divergence of whales. In ad- 
dition, we found that a model allowing a variation in co be- 
tween toothed whales and baleen whales is significantly 
better than a simpler model estimating the same co for both 
major lineages of whales (P=0.04, x 2 test). This result sug- 
gests differential levels of selective pressure acting on toothed 
and baleen whales. 

We attempted to identify all 24 coding exons of Trpm5 and 
all 32 coding exons of Plc/32 in the three whales with available 
genome sequences (fig. 5, supplementary table S2, 
Supplementary Material online). In the three Trpm5 se- 
quences, we detected multiple premature stop codons and/ 
or frame-shifting mutations, which are indicative of a Trpm5 
pseudogene in each of the three whales (supplementary table 
S2, Supplementary Material online). Notably, a 4-bp ORF-dis- 
ruptive deletion in exon 4 was shared between a toothed 
whale (Yangtze River dolphin) and a baleen whale (minke 
whale), suggesting the possibility of pseudogenization of 
Trpm5 in the common ancestor of whales. Additionally, we 
observed a 1 -bp ORF-disruptive deletion in exon 9 shared in 
the two toothed whales (bottlenose dolphin and Yangtze 
River dolphin). In contrast to Trpm5, Plcfi2 was found to be 
complete and intact in the three whales (supplementary table 
S2, Supplementary Material online), which indicated a func- 
tional role. To determine whether the relaxation of functional 
constraint occurred along the ancestral lineage of whales, we 
conducted a series of selection tests on Trpm5 and Plcfi2 that 
are similar to those on Calhml (table 1) using phylogenetic 
trees as shown in supplementary fig. S1C {TrpmS) and D 
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Calhml Exon 1 

Mus mus cuius (19 

Tursiops truncatus (19 

Stenella attenuate (19 

Lagenorhynchus albirostris (19 

Lagenorhynchus acutus (19 

Neophocaena phocaenoides (19 

Delphinapterus leucas (19 

Lipotes vexillifer (19 

Balaenoptera edeni (19 

Balaenoptera omurai (19 

Balaenoptera physalus (19 

Balaenoptera acutorostrata (19 

Calhml Exon 2 

Mws musculus (730 

Tursiops truncatus (730 

Stenella attenuate (730 

Lagenorhynchus albirostris (730 

Lagenorhynchus acutus (730 

Neophocaena phocaenoides (730 

Delphinapterus leucas (730 

Lipotes vexillifer (730 

Balaenoptera edeni (730 

Balaenoptera omurai (730 

Balaenoptera physalus (730 

Balaenoptera acutorostrata (730 



5 ' ATCTTCCAGTTCTTGCAATCCAACCAAGAGTCCTTCATGAATGGCATCTGTGGCATCATG 3 ' 

5 ' ATCTTCCAGTTCCTGCAGTCCAACCAGGAGTCCTTCATGAACGGCATCTGCAGCATCATG 3 ' 

5 ' ?????????????????????????????????????????CGGCATCTGCAGCATCATG 3 ' 

5 ' ??????????????????????????????TCCTTCATGAACGGCATCTGCAGCATCATG 3 ' 

5 ' ACCTTCCAGTTCCTGCAGTCCAACCAGGAGTCCTTCATGAACGGCATCTGCAGCATCATG 3 ' 

5 ' ATCTTCCAGTTCCTGCAGTCCAAGCAGGAGTCCTTCATGAACGGCACCTGCAGCATCATG 3 ' 

5 ' ATCTTCCAGTTCCTGCAGTCCAAGCAGGAGECClrTCATGAATGGCACCTGCGGCATCATG 3 ' 

5 ' ATCTTCCAGTTCCTGCAGTCCAACCAGGAGTCCTTCATGAACGGCATCTGCAGCATCATG 3 ' 

5 ' ATCTTCCGGTTCCTGCAGTCCAACCAGGAGCCCTTCATGAACGGCATCTGCAGCATCATG 3 ' 

5 ' ATCTTCCGGTTCCTGCAGTCCAACCAGGAGCCCTTCGTGAACGGCATCTGCAGCATCATG 3 ' 

5 ' ATCTTCCGGTTCCTGCAGTCCAACCAGGAGCCCTTCATGAACGGCACCTGCAGCATCATG 3 ' 

5 ' ATCTTCCGGTTCCTGCAGTCCAACCAGGAGCCCTTCATGAACGGCATCTGCAGCATCATG 3 ' 



5 ' GTATGTATCCAGCAGTTCTTTGAAGCCATG-AACCATGACCTGGAACTGGGTCATACCCAC 3 ' 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAC V 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAC 3 ' 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAT 3' 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAT 3' 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGTCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAT 3 ? 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAT 3 ' 

5 ' GTCTGTATCCAGCAGTTCTTCGAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAT V 

5 ' GTCTGTATCCAGCAGTTCTTCAAGGCCATG-ATCCATGACCTGGAGCTGGGTCACAGCCAT 3 ' 

5 ' gtctgtatccagcagttcttcaaggccatgEaSccatgacctggagctgggtcacggacat V 

5 ' GTCTGTATCCAGCAGTTCTTCAAGGCCATG^ATCCATGACCTGGAGCTGGGTCACGGCCAT 3' 

5 ' GTCTGTATCCAGCAGTTCTTCAAGGCCATG-ATCCATGACCTGGAGCTGGGTCACGGCCAT 3 ' 



Fig. 3. — An alignment of Calhml containing the ORF-disrupting mutations in seven toothed whales and four baleen whales. Dashes indicate alignment 
gaps, and question marks represent unamplif ied nucleotides. Codons in the correct reading frame are indicated by shading, whereas codons containing one 
nucleotide deletion or insertion are boxed. 



(Plcp2), Supplementary Material online. First, we analyzed 
data set III, which consisted of 14 nonwhale mammals and 
the ancestral sequence of all whales. We found that model F, 
which allowed a variation in co between the ancestral branch 
and other branches, is not significantly better than model E, 
which features a uniform co across the tree (P= 0.203). 
However, under the assumption of differential co ratios be- 
tween the ancestral branch leading to whales and the four 
branches connecting the three whales (model H), we found 
that model H to fit data set IV significantly better than model 
G, which assumes no co variation for these five branches 
(p= 0.005). Hence, these findings suggest that the relaxation 
of the functional constraint of Trpm5 did not originate in the 
ancestral branch of whales but may instead have occurred 
after the divergence of whales. Second, we examined data 
set V of 1 4 Plcfi2 sequences and found the ancestral branch of 
all whales to have a significantly higher co than other mammals 
(p= 0.007) after comparing model J with model I (table 1). 
Nevertheless, we did not find model L to be significantly better 
than model K (P=0.26), which indicates similar levels of se- 
lective pressure on the P/c/32 between the ancestral branch 
and branches connecting the three examined whales. 
Therefore, the relaxation of the functional constraint on 
P/c/32 may have occurred in the common ancestor of 
whales, although we failed to identify ORF-disruptive 



mutations in the sequences. In addition, we did not find sig- 
nificant differences in co between toothed and baleen whales 
for Trpm5 and PI02 (P> 0.4). 

Collectively, our evolutionary analyses suggest that the re- 
laxation of functional constraints on taste signal transduction 
likely took place in the common ancestor of whales between 
36 and 53 Ma, possibly shortly before 36 Ma, because the 
relaxed selection is weak, as indicated by small co ratios along 
the ancestral branch leading to whales (table 1). Thus, the 
relaxation of functional constraints on taste signal transduc- 
tion must have occurred after the pseudogenizations of taste 
receptor genes, probably because of other functions involved 
in these transduction pathways (Damak et al. 2006; Taruno 
et al. 2013). In addition, this relaxation is not a part of any 
genome-wide features, because the molecular evolutionary 
rate has generally slowed down in whales (McGowen et al. 
2012; Zhou et al. 2013; Yim et al. 2014). 

Discussion 

In this study, we surveyed all five primary taste receptor genes 
in toothed and baleen whales to examine the generality of 
whale taste loss and test whether taste loss has influenced the 
evolution of taste signal transduction. With the exception of 
the intact T2R16 in the three baleen whales, we found all 
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Tursiops truncates 
Stenella attenuata 
Lagenorhynchus albirostris 
Lagenorhynchus acutus 
Neophocaena phocaenoides 
Delphinapterus leucas 
Lipotes vexillifer 

91/1 00| Balaenoptera edeni 

Balaenoptera omurai 
Balaenoptera physalus 
Balaenoptera acutorostrata 



Myotis lucifugus 

Pteropus vampyrus 



Mustela putorius furo 




Ailuropoda melanoleuca 

Canis lupus familiaris 
Felis catus 
Oryctolagus cuniculus 



- Ochotona princeps 



Homo sapiens 
Pan troglodytes 
Gorilla gorilla gorilla 
Pongo abelii 
Nomascus leucogenys 



Macaca mulatta 

Callithrix jacchus 
Otolemur garnettii 



- Cavia porcellus 



Ictidomys tridecemlineatus 

Dipodomys ordii 

100/1001 Mms musculus 

Rattus norvegicus 



Loxodonta africana 



Fig. 4. — The ML tree of Calhml in mammals under the GTR+I+G substitution model of sequence evolution. Branch lengths are drawn to scale. ML 
bootstrap values/Bayesian posterior probabilities (>50%) are shown as numbers above the branches. 



members of sweet/umami (3 T1Rs) and bitter (1 0 T2Rs) recep- 
tor gene families (if amplified) to be pseudogenes in all 12 
examined whales, with a shared premature stop codon in 10 
of 13 genes, demonstrating massive losses of taste receptor 
genes in the common ancestor of toothed and baleen whales. 
Furthermore, we found that the sour taste marker gene 
Pkd2l1 is a pseudogene, whereas the candidate salty taste 
receptor genes {Scnnla, Scnnlb, and Scnnlg) are intact and 
evolutionarily conserved in two toothed whales and one 
baleen whale with available genome sequences. Finally, we 
showed the relaxation of selective constraints on taste signal- 
ing pathways in the common ancestor of whales (Calhml and 
Plc/32) or after the divergence of whales (Trpm5). Our findings 
unambiguously suggest widespread losses of umami, sweet, 
bitter, and sour tastes and the evolutionary conservation of 
salty taste in whales. The major reduction of taste receptor 
genes in whales may have resulted in the relaxation of selec- 
tive constraints on taste signaling pathways. Earlier studies 
discovered the loss of one or two primary tastes in certain 
groups of animals; our genetic study suggests that whales 



represent the first group of animals that strikingly lack four 
of five primary tastes. 

Consistent with our genetic evidence, anatomical studies 
have revealed taste buds in small pits at the base of dol- 
phin tongues and few taste bud-like structures in the ca- 
nonical taste structures (lingual papillae) of the four 
toothed whales (Yoshimura and Kobayashi 1997). 
Furthermore, taste buds were not discovered in seven 
other toothed whales (Arvy and Pilleri 1970; Yamasaki 
et al. 1976; Kuznetzov 1990; Pfeiffer et al. 2001). 
Behavioral tests of whale taste sensation are scant, with 
a few studies having examined bottlenose dolphins 
(Nachtigall and Hall 1984; Friedl et al. 1990). The dolphins 
could detect sweet, bitter, sour, and salty tastants with an 
order of magnitude below human sensitivities (Friedl et al. 
1990) or show indifference or reduced sensitivity to sweet 
and bitter tastants (Kuznetsov 1974). These anatomical and 
behavioral examinations could have been biased because of 
small sample sizes and few comparisons of sex and age 
differences. However, the morphological characteristics of 
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Table 1 

Likelihood Ratio Tests of Selective Pressures on Calhml, Trpm5, and Plcp2 in Mammals 



Models 


» (d N /d s ) 




Comparisons 


P 


Data set 1: 27 Calhml sequences (26 nonwhale mammals plus the ancestral sequence of all whales) 










A. All branches have the same co 


co = 0.06 








B. Ancestral branch of all whales has co 2 and other branches have cc^ 


co-i = 0.06, co 2 = 


0.17 


B vs. A 


0.021 


Data set II: 37 Calhml sequences (26 nonwhale mammals plus 11 whales) 










C. Ancestral branch of all whales and branches connecting 11 whales have co 2 , whereas other 


co-i = 0.06, co 2 = 


0.60 






branches have cot 










D. Ancestral branch of all whales has co 3i branches connecting 1 1 whales have co 2 , and other 


coi = 0.06, co 2 = 


0.69, 


D vs. C 


0.010 


branches have cot 


©3 = 0.16 








Data set III: 15 Trpm5 sequences (14 nonwhale mammals plus the ancestral sequence of all whales) 










E. All branches have the same co 


co = 0.09 








F. Ancestral branch of all whales has co 2 , and other branches have co-i 


(0-, = 0.09, (D 2 = 


0.15 


F vs. E 


0.203 


Data set IV: 17 Trpm5 sequences (14 nonwhale mammals plus 3 whales) 










G. Ancestral branch of all whales and branches connecting 3 whales have co 2i whereas other 


0^ = 0.09, co 2 = 


0.37 






branches have coi 










H. Ancestral branch of all whales has co 3i branches connecting three whales have co 2 , and other 


COt = 0.09, (0 2 = 


0.44, 


H vs. G 


0.005 


branches have cot 


co 3 = 0.12 








Data set V: 14 Plcft2 sequences (13 nonwhale mammals plus the ancestral sequence of all whales) 










I. All branches have the same co 


co = 0.11 








J. Ancestral branch of all whales has co 2 , and other branches have co-i 


coi = 0.1 1, co 2 = 


0.25 


J vs. I 


0.007 


Data set VI: 16 Plcfi2 sequences (13 nonwhale mammals plus 3 whales) 










K. Ancestral branch of all whales and branches connecting three whales have co 2 , whereas 


coi = 0.1 1, co 2 = 


0.19 






other branches have coi 










L. Ancestral branch of all whales has co 3; branches connecting three whales have co 2 , and 


co-i = 0.12, co 2 = 


0.16, 


L vs. K 


0.254 


other branches have 


co 3 = 0.25 









Note. — Significant P values (< 0.05) are indicated in bold. 



Trpm5 Exon 1 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(58) 5' GAGCCCATCCTATGCAGGGGAGAGATCAACTTCGGAGGGTCTGGGAAGAAGCGAGGCAAG 3' 

(58) 5' GGGCCGGGCCTGTGCCAGGGTGAGATCGAATTTGGAGGGTCTGGGAAGAAGCGAGGCCAG 3' 

(58) 5' GGGCCGGGCCTGTGCCGGGCTGAGATCGGATTCGGAGGGTCTGGGAAGAAGlGlGGCCAG 3' 

(58) 5' GGGCCGGGCCTGTGCCGGGGTGAGATCGAATTCGGAGGGTCCGGGAAGAAGCGAGGCAAG 3' 



Trpm5 Exon 4 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(559) 
(559) 



CTTGTGGAGTCAGGCGCCCTTGGGAGTGGGAACGACGGGCTGACAGAGCTGCAGCTGAGC 

999999999999999999999999999999999999999999999999999999999999 



(559) 5' CCGGTGGACCCGGGCACCCCTGGGAAGEE3E3CACCCGCCCGCAGAGCTGCGGCCA 

(559) 5 ' CTGGTGGACCCGGGCGCCCTGGGGAGGE3E3AGCGAGCCCACTGAGCTGCGGCTGAGG 



Trpm5 Exon 8 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(1078) 5 ? GCCTGGGATCGCGTGGACATTGCCAAGAGTGAAATCTTCAATGGGGACGTGGAATGGAAG 3' 

(1078) 5' GCCTGGGACCGCGTGGACATCGCCAGGAGCGAGATCOAACAGGGACGTGCAGTGGAAG 3' 

(1078) 5' GCCTGGGACCACGTGGACATCACCAGGAGCGAGATCTTCAACGGGCACGTGCAGCGGAAG 3' 

(1078) 5' GCCTGGGACTGCGTGGACATCGCCAGGAGCGAGATCCTCAACGGGGACGTGGAGTGGAAG 3' 



Trpm5 Exon 9 
Mouse 

Bottlenose dolphin 
Yangtze River dolphin 
Minke whale 



(1207) 5' gtggacagcggtgctgacatggccgagttcttgacctatgggcggctgcagcagctttac 3' 

(1207) 5' gtggacagtggcgccaacgtggcggacttcEi3acgtacgggcggctgcagcagctctac 3' 

(1207) 5' gtggacagtggcgccgacgtggcggacttcFtgIacgtacgggcggctgcagcagctctac 3' 

(1207) 5' gtggacggtggcgccgacgtggcggacgtcgtgacgtacgggcggctgcggcagctctac 3' 



Fig. 5. — An alignment of Trpm5 containing the ORF-disrupting mutations in two toothed whales (bottlenose dolphin and Yangtze River dolphin), one 
baleen whale (mink whale), and the outgroup (mouse). 
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the whale tongues are generally simplified, and their taste 
sensitivity is poorly developed. 

Why could whales afford to lose four of five primary tastes? 
We provide three probable explanations. First, the feeding 
behavior of swallowing prey whole without mastication may 
have rendered their tastes useless, as suggested by two recent 
studies (Jiang et al. 2012; Sato and Wolsan 2012). Second, 
their food items are likely to possess reduced taste stimuli 
because a high concentration of sodium in the ocean could 
have masked the tastant cues, which could be quickly diluted 
in sea water (Ikeda 1909; Komata 1990). Third, a dietary 
switch from plants to meat in the whale ancestor may account 
for the major loss of sweet and bitter tastes, because meat 
contains little sweet and bitter compounds. This scenario is 
analogous to that occurred in the giant panda, which has 
lost its umami taste due to a dietary change from meat to 
bamboo (Zhao, Yang, et al. 2010). The extinct even-toed un- 
gulates (Raoellidae) that are most closely related to whales 
were plant eaters, which strongly suggests that a major dietary 
change occurred during the transition from even-toed ungu- 
lates to whales (Thewissen et al. 2007). 

In contrast, the maintenance of the only primary taste, salty 
taste, is not unexpected. The sodium channel ENaC, which 
functions as the candidate salty taste receptor, typically con- 
sists of a multimeric complex of three subunits and plays a 
major role in the regulation of extracellular volume and blood 
pressure (Canessa et al. 1994). Each of the three subunits 
showed significant sequence similarity and similar expression 
patterns between invertebrates and mammals (Chang et al. 
1996; Ottaviani et al. 2002). This evolutionary conservation of 
ENaC strongly suggests its functional indispensability across 
animals. Indeed, the hereditary mutations of ENaC subunits 
were demonstrated to result in Liddle's syndrome of hyper- 
tension (Schild et al. 1996); disruptive mutations creating 
dysfunction of ENaC were found to cause pseudohypoaldos- 
teronism with severe sodium loss (Chang et al. 1996; Kerem 
et al. 1999). As sea-living organisms, whales are well adapted 
to their hyperosmotic environment; osmoregulation is re- 
quired to maintain the ENaC function that regulates sodium 
ion fluxes and thus retains the salty taste (Ottaviani et al. 
2002). 

Taste receptor evolution in relation to feeding ecology is 
not well understood. The evolutionary patterns of taste recep- 
tor genes matched the variations of feeding ecology in the 
giant panda (Zhao, Yang, et al. 2010), humans (Wang et al. 
2004), vampire bats (Zhao et al. 2010), carnivores (Jiang et al. 
2012), and vertebrates in general (Li and Zhang 2014), but 
mismatches between taste receptor evolution and feeding 
ecology were also observed (Zhao and Zhang 2012; Zhao 
et al. 2012). Thorough scrutiny of the physiological functions 
of taste receptor genes should help to understand the con- 
nection between taste function and feeding ecology. Indeed, 
in addition to a gustatory system, taste receptors are also in- 
volved in extragustatory systems, such as the lung, gut, brain, 



and testis (Behrens and Meyerhof 201 1), which suggests that 
these extragustatory functions may account for the aforemen- 
tioned mismatches. Nonetheless, the massive losses of taste 
receptor genes in whales are striking, although losses of taste 
receptors have been frequently identified from insects to 
mammals (Nei et al. 2008; Xiao et al. 2013). Therefore, 
taste receptors that are pseudogenized could not be involved 
in extraoral functions in whales; the functional significance of 
these extraoral taste receptors in humans and rodents must be 
unimportant in these marine mammals. Future studies of 
other sensory systems that are still maintained in whales 
may provide a better understanding of how whales and 
other marine animals can sense and survive in the ocean. 

Supplementary Material 

Supplementary figure S1, tables S1 and S2, and data set S1 
are available at Genome Biology and Evolution online (http:// 
www.gbe.oxfordjournals.org/). 
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